Solution of the local field equations for self-generated glasses 
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We present a self-consistent local approach to self generated glassiness which is based on the 
concept of the dynamical mean field theory to many body systems. Using a replica approach to 
self generated glassiness, we map the problem onto an effective local problem which can be solved 
exactly. Applying the approach to the Brazovskii-model, relevant to a large class of systems with 
frustrated micro-phase separation, we are able to solve the self-consistent local theory without using 
additional approximations. We demonstrate that a glassy state found earlier in this model is generic 
and does not arise from the use of perturbative approximations. In addition we demonstrate that 
the glassy state depends strongly on the strength of the frustrated phase separation in that model. 

PACS numbers: 



I. INTRODUCTION 

The familiar example of vitrification upon super- 
cooling molecular fluids provides proof that an effectively 
non-ergodic disordered state can be generated in a sys- 
tem without pre-existing quenched disorder. Such self- 
generated glassiness may be a more widespread possibil- 
ity in the material world than is currently acknowledged. 
Candidate solid state systems with extremely slow dy- 
namics often have impurities remaining from their syn- 
thesis so there is always a suspicion that non-ergodicity 
comes from di sor cleri2i2i2i£i§iL£i2ii2iiL While a rather 
successful microscopic theory of molecular glasses now 
exists based on the idea of an underlying random first or- 
der transitioni£ii§*iLi&, this theory, like all theories of liq- 
uid state phenomena, contains approximations that can 
only be checked a posteriori through their agreement 
with simulations or experiments. To deepen and widen 
our understanding of glassy dynamics beyond the molec- 
ular fluid example, it would be desirable to have a model 
which can be analyzed in a formally exact fashion using 
familiar tools of continuum field theory. Such a theory 
would also help clarifying which features of the current 
theory of structural glasses are robust and "protected"— 
and which are fortuitously correct. 

Both, the recent treatment of the thermodynamics of 
fragile glasses developed within a replica formulation 13 
and the closely related approach using density functional 
theory for super-cooled liquids^ suggest that the major 
aspects of glass formation stem from the very strong 
local correlations of a dense super-cooled fluid. This 
is very different from systems with critical dynamics 
close to second order phase transitions, where long wave- 
length fluctuations are paramount. One finds experimen- 
tally that dynamics on essentially all wavelengths slows 
to nearly the same extent as the glass transition is ap- 
proached. This suggests that the search for an solvable 
limit of the glass problem should focus on the strong but 
local nonlinearities. 



Major progress has recently been made in the the- 
oretical understanding of strongly interacting systems 
with predominantly local correlations. Strongly corre- 
lated electron systems like Mott insulators or systems 
with strong electron-phonon interactions have been inves- 
tigated with the dynamic mean field theory 36-37 (DMFT). 
Originally, this approach was based on the assumption 
that all non-trivial correlations of a given system are 
strictly local (for the precise definition of a local theory, 
see below) and then allowed for a solution of such sys- 
tems without further approximations. More recently sev- 
eral generalizations of the approach have been developed, 
which allow the inclusion of collective excitations^^ as 
well as short range correlations in so called cluster DMFT 
theoriea 23 : 24 : 25 . 

In this paper we will follow the main strategy of the 
dynamic mean field approach and apply it to the replica 
theory of self generated glassiness. This allows us to 
solve completely for the long time correlation function 
(Edwards- Anderson parameter) as well as the configu- 
rational entropy of a glass within a local theory. Even 
though the origin of glassiness is not quenched disorder, 
there is a similarity between our theory and the mean 
field theory of spin glasses which is also exact in the 
limit of a purely local approach. Solvable local limits 
like the Shcrrington-Kirkpatrick model, turned out to be 
extremely fruitful for the investigation of more realistic 
models. To identify and solve such a local problem is 
therefore of general interest. Another appealing aspect of 
the approach developed here is that the recent generaliza- 
tions of the DMFT to cluster systems offer a feasible way 
to develop a controlled series of approximations for can- 
didate glassy systems that successively extend the range 
of local correlations taken into account— Clearly, 
before this becomes possible we have to set up and an- 
alyze the strictly local (i.e. mean field) theory. The so- 
lutions of the local field theories discussed in this pa- 
per rely on a self consistency condition that is uniformly 
applied to the replicated statistical mechanics problem. 
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The droplet effects that give rise to the activated tran- 
sitions in molecular fluids cannot be described by this 
uniform replica symmetry breaking ansatz. While the 
self consistent replica field theory is applicable to account 
for some geometric effects in finite dimensions the latter 
highly nonperturbative effects are absent, but are sup- 
pressed if the coordination number is high. We should 
also mention that related interesting approach to self- 
generated glassiness in a model in infinite dimensions 
(where the locality of the problem is exactly guaranteed) 
that was recently presented by Lopatin and IofTe^ 

We demonstrate that a uniformly frustrated system 
with competition of interactions on different length scales 
undergoes a glass transition once inhomogeneous fluctua- 
tions become strong enough if no ordered phase forms via 
nucleation. Specific calculations are carried out for the 
case of the Brazovskii model 4 -, which describes systems 
with tendency towards micro-phase separation. Glassi- 
ness was shown to exist in this model in Refs l2al3ll How- 
ever, it remained unclear to what extent this is merely a 
consequence of an approximate solution using perturba- 
tion theory or is indeed the correct mean field behavior 
of this model. 

In the next section we briefly summarize the replica 
approach to self generated glassiness, present the main 
idea of the dynamic mean field theory and solve the local 
theory for self generated glassiness in case of the Bra- 
zovskii model. We summarize the results of this paper 
in Section III. 



II. THEORY 

A. The replica approach 

In this paper we consider physical systems consisting 
of L-components (types of atoms or molecules) and with 
short ranged 3 rd , 4 th and higher order virial coefficients. 
No specific assumption is made with respect to the (pos- 
sibly non local) 2 nd virial coefficient. In terms of a 
density field pi (x) = £\ <5 (x — X^;) of particles of type 
I = 1, L ( L = 1 hx case of a one component system) 
we then have an effective energy: 

H[4>]= H [</>]+ [ d d xV (1) 



with 

H = [ d d xd d x' 4*1 ( x ) G^], (x, x') fa, (x') . (2) 
^ w 

where <pi (x) = pi (x) — ~p l is the deviation of the den- 
sity from its average value. The key assumption made 
in EqsQ] and [21 is the local character of V (<p) which in- 
cludes the higher order virial coefficients. For example, 
Eq^ might be a representation of the density functional 



used by Ramakrishnan and Youssoff with Gq W (x, x') = 
cui (x — x') the direct correlation function of the fluid 35 
and 

v (4>) = fa + h (pi + <h (x)) - 1] (3) 

the ideal gas free energy. 

The partition function of the system is given by 

Z cq = J ityexp(-ff M/T) (4) 

and determines the equilibrium properties of the system. 
In a system with complex energy landscape where we 
expect that the system might undergo vitrification, the 
knowledge of the equilibrium partition function is not 
sufficient anymore. 

A widely accepted viewi***^ is that a glassy system 
may be considered to be trapped in local metastable 
states for very long time and can therefore not realize a 
considerable part of the entropy of the system, called 
the configurational entropy S c — log J\f m a, where J\f m s 
is the number of metastable states. If Af m s is exponen- 
tially large with respect to the size of the system, S c 
becomes extensive and ordinary equilibrium thermody- 
namics fails. There are several theoretical approaches 
which offer a solution to this breakdown of equilibrium 
many body theory. On the one hand one can solve for the 
time evolution of correlation and response functions, an 
approach which explicitly reflects the dynamic character 
of the glassy state. Mostly because of its technical sim- 
plicity, an alternative (but equivalent) approach is based 
on a replica theorjiiSil 3 .. Even though this approach docs 
not allow one to calculate the complete time evolution, 
long time correlations as well as stationary response func- 
tions can be determined which are in agreement with the 
explicit dynamic theory. We will use the replica approach 
because of its relative simplicity. 

The central quantity of the replica theory to self gen- 
erated glassiness is the replicated partition functioniSii 3 . 

Z (m) = J D ro ^sr ff iw +s/^e:>°( x )m x ). 

(5) 

Z (m) must be analyzed in the limit m —> 1 and g — + 
which is taken at the end of the calculation. One can in- 
terpret Z (to) as the result of a quenched average over 
random field configurations with infinitesimal disorder 
strength g and a distribution function which is deter- 
mined by the partition function of the system itself. A 
similar approach (however for finite g) was proposed by 
Deam and Edwards^ in the theory of the vulcanization 
transition where the distribution of random cross-links of 
a polymer melt is also determined by the actual partition 
function of the polymer system itselfiS. Then, just as in 
the present approach, the number of replicas, to, is taken 
to the limit m — > 1, as opposed to the usual m — > limit 
used for systems with averaging over "white" quenched 
disorder distributions. In the vulcanization problem, g 
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is related to the cross link density and, in distinction to 
the theory presented here, remains finite. Physically, the 
analysis of Z (m) addresses whether the system drives it- 
self (m — ► 1) into a glassy state if exposed to an infinites- 
imally weak randomness (g — ► 0). It can be shown^ that 
this approach is equivalent to the solution of the dynamic 
equations for the correlation and response function^ of 
the system that takes aging and long term memory effects 
into account. 

F (m) = — ^ log Z (m) determines the configurational 
entropy, S c , and the averaged free energy of systems 
trapped in metastable states, F, viai2*i^i 



F 



lim 
lim 



d (mF (m)) 



dm 
i 2 OF (m) 



T dm 



(6) 



The equilibrium free energy, F = — T log Z cq can be 

shown to be F = F — TS C - This makes evident that 
TS C is indeed a part of the free energy which cannot be 
realized if the system is trapped in metastable states. 

The limit m — > 1 is only appropriate above the Kauz- 
mann temperature, Tk, at which S c vanishes. Below Tk 
one has to perform the limit m — > T/T e g with effective 
temperature T e g > T characterizing the spectrum of 
metastable states. How to determine T e g in case of a 
rapid quench is discussed in Ref. EE whereas an inter- 
esting alternative to determine T e g for slowly quenched 
systems was proposed in Refill In the present paper we 
restrict ourselves to the behavior T > Tk even though 
our approach can be generalized to the case with arbi- 
trary m < 1 easily. 



B. Mapping onto a local problem 

The major problem is the determination of the parti- 
tion function Z (m). Even for the liquid state (i.e. m = 1 
and g = at the outset) this is a very hard problem 
without known exact solution and we are forced to use 
computer simulations or to develop approximate analyt- 
ical theories. In developing such an approximate theory 
we take advantage of the fact that glass forming systems 
are often driven by strong local correlations, as opposed 
to the pronounced long ranged correlations at a second 
order phase transition or the critical point of the liquid- 
vapor coexistence curve. This is transparent both, in 
the mode coupling theory of under-cooled liquids 34 and 
in the self-consistent phonon approaches 1 *, where a given 
molecule is locally caged by its environment built of other 
molecules. 

When we call a physical system local this does not nec- 
essarily imply that all its correlation functions infinitely 
rapidly decay in space. In the language of many body 
theory it only implies that the irreducible self energy 
(k) ~ Y>ui is independent of momentum for the im- 
portant range of parameters. Here, £;;' is related to the 



correlation function Gu> (k) 
equation 



fy,k#,-k) via Dyson 



G(k)" 1 =G (kV^-E 



(7) 



which is an L x L matrix in case of the equilibrium liq- 
uid state theory. If we study the emergence of glassy 
states we have to use the replica theory and Eq0 be- 
comes an (mL x mL) matrix equation with G^f (k) = 

Ot'l^V ,-k)> G ou> ( k ) = G ow (k)<W as well as self en- 
ergy matrix E^f. 

Traditionally, the self energy is introduced because it 
has a comparatively simple structure within perturba- 
tion theory. However, in the theory of strongly correlated 
Fermi systems it has been recognized that the existence of 
a momentum independent self energy allows conceptually 
new insight into the dynamics of many body systems, de- 
scribing states far from the strict perturbative limit 
We will adopt the main strategy of this dynamic mean 
field theory for our problem. 

We use the fact that the free energy F (m) determined 
by EqJSJcan be written as2£ 

F (m) = Itr log (GoT 1 - S) + hr (EG) + [G] , (8) 

where the trace goes for each k-point over the 
(mL x mL)-matrix components together with a sum over 
k. The latter can also be written as a matrix trace of 
real space functions G(x, x') etc. The functional $ [G] 
is well defined in terms of Feynman diagrams as the sum 
of skeleton diagrams of the free energy. In what follows 
we will not try to calculate $ but merely use the fact 
that such a functional exists. From the definition of $ it 
follows that it determines the self energy via a functional 
derivative: 



E«f(x,x') 



S<S> [G] 



SGft (x',x) 



(9) 



Since we made the assumption that £ is momentum 
dependent, Fourier transformation yields E^f (x, x') = 
£^f <5 (x — x'). Thus, the functional derivative, Eq03 
vanishes if x ^ x', which implies that for a local the- 
ory, <I> solely depends on the local, momentum averaged, 
correlation function, 



(10) 



Since all our interactions, V [4>] are (by assumption) local 
as well, we conclude that there exists a local problem with 
Hamiltonian 

ll',a/3 a 

which has an identical functional <I> [Q] of its own cor- 
relation function Q, which is also a [mL x m_L)-matrix 
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but does not depend on position or momentum, ao is a 
typical microscopic length scale, for example a hard core 
diameter and needs to be specified for each system. It 
results from the fact that a sensible local theory can only 
be formed after some appropriate discretization 



d d xV(fr{x))^a$J2 V (Mxi))- (12) 



Even though, H. has no spatial structure anymore, the 
perturbation theory up to arbitrary order is the same for 
both systems. This holds for an arbitrary choice of the 
so called Weiss field, J. 

Following Refl37l we use the freedom to chose J in 
order to guarantee that Q = G. This implies that not 
only the functional $ but also its argument are the same 
for the actual physical system and the auxiliary local 
one. It then follows that the self energy of the original 
system E^jf is, up to a trivial prefactor which results 
from the above discretization, equal to the self energy of 
the auxiliary system, a^E^f . The prefactor ao can for 
example be determined by comparing the leading order 
in perturbation theory of the two problems, EqQ and 
Ea llll Then we solely need to solve the much simpler 
problem, TL 1 and determine for an assumed Weiss field 
J the local self energy as well as the local propagator 
related by 



G = (J - ajfE) 



(13) 



We made the right choice for J if simultaneously to EalTSI 
it is true that self energy and averaged correlation func- 
tion are related by: 



(14) 



If this second equation is not fulfilled we need to improve 
the Weiss field J until Eg Upland Eq^Jhold simultane- 
ously, posing a self consistency problem. This is the most 
consistent way to determine the physical correlation func- 
tions under the assumption of a momentum independent 
self energy. This approach has been applied to a large 
class of problems in the field of strongly correlated Fermi 
systems, and promises, as we argue here, to be very use- 
ful in rather different contexts. The major technical task 
of the DMFT is to solve H for given J. This will be done 
for the specific choice of a replica symmetric correlation 
function 



Gtf (q) - Kui (q) 5 a0 + Fui (q) , 
which implies a similar form for the self energy 
E Q/3 = Y, HI 5 a p + 



(15) 



(16) 



Below we will present a general approach to study the 
stability of this choice and analyze it for a specific ex- 
ample. 



C. Example: the Brazovskii model of micro-phase 
separation 

We are now in the position to apply our approach to 
a specific physical system. We consider a one compo- 
nent system (L = 1) governed by the three dimensional 
Brazovskii model^, 



H=Ud*X Uqlcf + \<f 



(17) 

which has a broad range of applicability in systems 
with micro-phase separation like the theory of micro- 
cmulsiona^2iMi2&, block copolymers " or e ven doped 
transition metal oxidesS&^iili. In Refs l40l4ll it was ar- 
gued that it might be used as a simple continuum model 
for glass forming liquids. The bare correlation function 
follows from EcillTI 



Go (?) 



1 



4ll 



1 



2 2 

4% 



(q - qoY 



(18) 



The interesting behavior of the Brazovskii model arises 
from the large phase space of low energy fluctuations 
which is evident from the gradient term, V 2 + g 2 , in the 
Hamiltonian. All fluctuations with momenta |q| = go 
can, independent of the direction of q, be excited most 
easily. The wave number go is related to various physical 
quantities in all these different systems. In microemul- 
sions, go is determined by the volume fraction of am- 
phiphilic molecules whereas it is inversely proportional 
to the radius of gyration in block copolymers and to 
the strength of the Coulomb interaction in doped tran- 
sition metal oxides. Clearly, the role played by the 
microscopic length scale ao of the previous section is 
q^ 1 ■ As discussed in Appendix A, an explicit calcu- 
lation within the self consistent screening approximation 
gives ao — l.TqQ 1 . For simplicity we use ao = g^ 1 , which 
corresponds to a rescaling of the coupling constant. 

If the self energy is momentum independent the corre- 
lation function has the form 



G(q) = 



1 



(q - qo) 



e 2 g 2 ' 



with e 2 = e\ — q Q 2 E > 0, which yields 



G 



— TG(g) = — . 



(2tt 



(19) 



(20) 



In Refs l2i8ll it was shown that a self generated glass 
transition can occur and the configurational entropy was 
calculated. A glassy state occurs for all finite go and 
leads to S c oc VqQ, with volume V. However, these re- 
sults were obtained by using perturbation theory and in 
particular the discussion of Ref|29j suggests that the self 
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consistent screening approximation might strongly over 
estimate the tendency towards glassiness. An important 
ques tion is therefore whether the glassy state obtained in 
Refs l28l3ll is generic or arises solely as consequence of a 
poor approximation. This issue was also raised in the 
recent numerical simulations in The approach de- 

veloped in this paper is able to answer this question. In 
order to facilitate the comparison with the pertu rbative 
results we summarize the results of Refs. I28l3ll in Ap- 
pendix A however using an argumentation more adapted 
to the dynamic mean field approach. 



1. DMFT for the liquid state 

Ignoring glassiness for the moment, the local Hamilto- 
nian is given by 



(21) 



which leads to the partition sum Z = dipe n and 
correlation function 



G 



00 AAA2 e -U 



(22) 



Z and G are elementary integrals that can be expressed 
in terms of elliptic functions. The Weiss field which 
leads to the correct propagator can be obtained as J = 

e ° +2 qo~ £2 wnich > together with Eqs[!)| and EU leads to 
a single nonlinear algebraic equation for e. The solution 
can be obtained easily. Results for Eq = —4 (strong seg- 
regation limit) as well as £q = —1 (weaker segregation 
limit) are shown in Fig.l and are compared with the cor- 
responding behavior as obtained using a Hartree theory 
which gives EHartrcc = — -i^r ■ Leading order pertur- 
bation theory (i.e. Hartree theory) does give the proper 
small u behavior, however strong deviations occur in case 
of the weak segregation limit (small £q ) ■ 

In our treatment of the equilibrium behavior we have 
made the assumption that no phase transition to a state 
with long range order of (</> (x)) takes place. This is con- 
sistent with the Hartree theory which excludes a diverg- 
ing correlation length for finite temperatures, and thus, 
excludes the possibility of a second order phase transi- 
tion. However this assertion contrasts with the result 
found by Brazovski: 44 who showed that the model, Ea ll7l 
undergoes a fluctuation induced first order transition to 
a modulated state with order parameter 



> (x)) = (j) cos (q 



(23) 



where |qo| = qo and arbitrary but fixed direction. For 
small coupling constant this transition occurs if Eg > 



(i) 



2/3 



— - i — , The criterion for the glass transition 

(Sir) 2 '- 1 \1o J ° 

discussed in this paper is very similar to the parame- 
ter values leading to the fluctuation induced first order 



GO 



DMFT 

Hartree approximation 




FIG. 1: Dimensionless inverse correlation length of the equi- 
librium liquid state as function of the coupling constant 
(scaled by the modulation wave number qo for different segre- 
gation strength £q. The DMFT results (solid lines) are com- 
pared with the leading order Hartree approximation (dashed 
lines). Strong deviations for small dimensionless coupling con- 
stants it/qo occur in the limit of weak segregation strength. 
Note, the results are obtained by ignoring the possibility of 
a fluctuation induced first order transition which will become 
relevant if | e§ | > c (u/9o)'2/3) with c a constant of order unity. 



transition. We have not analyzed both transitions within 
an identical theoretical approach but find in Appendix A 
that perturbation theory yields for the occurrence of the 

3/5 



glass transition |sq| cx (jjj^ ■ bi what follows we will 

always assume that the glass transition occurs instead of 
the fluctuation induced first order transition. In the lab- 
oratory this will likely be a kinetic issue, that needs 
a nucleation theory for quantitative predictions. The 
results of Reffjjl where a nucleation theory of the Bra- 
zovskii transition was developed demonstrate that indeed 
the nucleation kinetics of the order parameter, Eg 1231 is 
very complex, i.e. under-cooling should be possible. 



2. DMFT in the glassy state 

To use the replica theory of the glassy state we have to 
specify the replica structure of the correlation function. 
We first choose a given structure and discuss its stabil- 
ity later. In replica space we start from the following 
structure of the propagators and self energies: 



G a p{q) = K (q) S al3 + F (q) 
S a /3 = TiKO-aP + Sf. 



(24) 
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Inverting the matrix Dyson equation in replica space 
leads to: 



K{q) 



mF (q) 



(q - Qo) + K 2 ql 
1 



(25) 



(q-Qo) + e 2 <7o (q-qo) + ^ 2 ql 

with the convenient definitions: 



2 2 

K q 

2 2 



-2„2 



(q) 



(26) 



e q - H K (q) ~ mT, F (q) . 



The diagonal elements K (q) + F (q) can be interpreted 
as the equilibrium, liquid state correlation function and 
is only determined by e which can be related to the liq- 
uid state correlation length f ~ — . On the other hand 
F (q) = linij^oo lim^oo (<fi q (t) (f>- q (t + t')) character- 
izes long time correlations in analogy to the Edwards- 
Anderson parameter. Clearly, if E^ ^ then k > e 
and F (q) > 0. Depending of whether one considers q 
values close to or away from qo, F (q) is governed by the 
correlation length £ or the Lindemann length of the glass 

Ao = q$ (n 2 — s 2 ) , respectively. The physical signif- 
icance of Ao is the length scale over which defects and 
imperfections of an crystalline state can wander after 
long time was discussed in Ref. 31 . Finally, the correla- 
tion function K (q) (which is solely determined by the 
short length < £) jg the response function of a local 
perturbation. Obviously, any response of the glassy sys- 
tem is confined to very small length scales even though 
the instantaneous correlation length can be considerable. 
This is a clear reflection of the violation of the fluctuation 
dissipation relation within the replica approach. Averag- 
ing these functions over momenta gives: 



K 
F 



^F{q) = ^(--~ 
(2tt) 3 2ttto \e k 



(27) 



The auxiliary local Hamiltonian is 



«=5£^a& + ^5>*. (28) 
ab a 

Within DMFT we then find for its correlation function 

g = K5 ab + F. (29) 

In addition, the Weiss field is given by: 

Jafl = J5 a/3 - C. (30) 

These relations together with EqE|can be used to ex- 
press the Weiss fields in terms of the e and k: 



J = 



2ttk + Eq — K 2 

qo 

2ir (k - e) + e 2 - k 2 
mq 



(31) 



Thus, we have to determine K and F for given J and C 
and make sure that the latter are chosen such that Easl2*7l 
and OH] are fulfilled. 

The partition sum of the local problem is given by 



z{m) = J «r^-E: =1 «o[^]+^ a ^) 2 



(32) 



refers to the fact 



where H [<j> a ] = \J4>\ + ' 
that 4> is an m-component vector and the integral goes 
over an m-dimensional space with arbitrary m. The cou- 
pling between different replicas can be eliminated by per- 
forming a Hubbard-Stratonovich transformation, which 
leads to 



where 



Z(m) 



ft (A) 



d<j)e 



-(Hal4>]+\<t>) 



(33) 



(34) 



is the equilibrium partition function however in an exter- 
nal field A, with Gaussian distribution function. 

In order to determine the propagators of the local prob- 
lem, we consider the sum of the diagonal elements 



^ G aa = m 



which is equal to —2 



d log Z{va) 
dj 



K + F) 



yielding 



(35) 



K + F = -2 



Z (m) V2nC dj 
The derivative with respect to J leads to 



n (A) . (36) 



d_ 

dj 



n (A) 



-U^ 2 +- 



-0 4 + A0 



4? (37) 



which gives the final expression for the diagonal element 
of the replica correlation function 



with 



K + F 



Z (m) 



(38) 



SI (A) 



(39) 



In addition we also need to determine the off diagonal 
elements in replica space of the correlation function. We 
use 



J2 G «b = mK + m 2 F 



(40) 
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which equals to 2 dl ° 6 g ^ m " > and obtain an equation which 
can be used to determine the off diagonal elements F 



K + mF 



dX 
/2nC 



Z(m) 



(41) 



Evidently, Eqs|35]andE]are only independent if m differs 
from 1. 

We next analyze these equations for small but finite 
m — 1. This can be done by expanding EaslSSI and FTTI 
into a Taylor series for small m — 1 and comparing order 
by order. First we consider the zeroth order term and 
find that both equations yield for m = 1, the results 
for the liquid state which determine e. The next step 
is to consider the first corrections linear in m — 1. This 
allows us to check whether there are nontrivial solutions 
for k > e and thus for the off diagonal self energy and 
long time correlation function. 

The difference between Easl38l and 1411 gives 



F 



(1 — m) Z (m) 
Expanding the numerator for m close to 1 gives 

SI 



F 



Z(m = 1) 



(42) 



(43) 



with 



51 = 



I 



d\ (C- 1 (4>) x A + <0 2 ) A ) e-^c O (A) log SI (A) 

(44) 

As expected, it follows that F = if n = e, where C = 
0. This can be seen from the expansion (0) A ~ — (</> 2 ) A, 

valid for small A and by substituting fi = A/VC. It fol- 
lows that 



F (k = e) cx 



d ^-_e-^ (l-/i 2 ) =0. (45) 



>2it 



In Fig. 2 we show our results for the dimensionless in- 
verse correlation length £a = £ (ma) at the onset of glassi- 
ness. We have to distinguish the two different regimes 
of strong and weak phase separation, characterized by 
|ej)| > 1 an d | e§ | < 1, respectively. In the case of strong 
phase separation we find a universal value £a — 0.45. 
The universality of this result is in qualitative agreement 
with results obtained with perturbation theory which 
also gives such a generic value which is however numer- 
ically larger (ea — 0.8385). In addition we obtain that 
at the transition the ratio ka/sa — 3.5, which is close 
to the result ka/ea — 3.1, obtained within perturbation 
theory. The onset of glassiness in the model EallTI occurs 
in the strong phase separation limit and, for given q , 
is uniquely determined by the value of the correlation 
length £ = (eqo) , and not explicitly dependent on u or 



0.7 



0.6 



0.5 



1 1 1 1 

liquid e A SCSA (u^oo)=o.84. 




FIG. 2: Dimensionless inverse correlation length as function 
of the coupling constant at the dynamical glass transition ea 
obtained within the DMFT, in comparison with the result 
of the self consistent screening approximation. For large u 
a universal value of the correlation length (in units of q§ 
for glassiness results, whereas the critical value of the spatial 
correlations to reach a glassy state increases considerably in 
the weak coupling limit. 



T . The exact solution of the dynamic mean field prob- 
lem leads to the interesting result that this value for £ 
is larger than what follows within perturbation theory. 
For small the result of the perturbation theory is 

e A a 0.42 If) .In this limit the DMFT gives the in- 
teresting result that £a is larger than the value obtained 
within perturbation theory. The effect is small and nu- 
merically hard to identify. It implies that the criterion 
for glassiness as obtained within DMFT becomes very 
close to the one found by Brazovskii for the fluctuation 
induced first order transition. Most important however 
is that we do indeed find that there is a nontrivial so- 
lution k > £. The emergence of a self generated glassy 
state is therfore not a consequence of the perturbative so- 
lution used in RefsHUll Our calculation is performed 
for a temperature T = 1. We can reintroduce the tem- 
perature into the calculation by substituting u — > uT . A 
critical value for u leads to a temperature Ta (for fixed u) 
where within mean field theory an exponential number 
of metastable states emerges. 

Once the correlation functions are determined we can 
use the fact that the functional $ [G] is the same for the 
local problem as well as the original one to obtain 

F(m)=J r (m) + —( tr logG" 1 - a^Tr logG" 
2m V 

(46) 

Here Tr refers to the trace over replicas, but does not in- 
clude the momentum integration, as opposed to tr which 
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D. Stability of the solution 




FIG. 3: Configurational entropy density S c /V (with system 
volume V) as function of temperature for two different values 
of the segregation strength 6q. 



corresponds to a trace with respect to all degrees of free- 
dom. Finally, T (m) = — ^ log Z (m) is the counter part 
of F (m) for the local problem. Using Eq|S|an expression 
for the configurational entropy follows: 



F(q)_ ) + F(q). 



.-d 



-(1-f) 



G(q)' G(q) 
_ d d log Z (m) /m 



dm 



(47) 



where G(q) = K(q) + F(q) and 



dZ (m) 



dm 



(A) log ft (A). (48) 



m — *1 



As expected, in the limit n — > £ without glassy long 
time correlations, it follows 2(1) = f2 (0) as well as 

9z ( m } = Q (0) log ft (0) and F(q) goes to zero, lead- 

m m— »1 

ing to (k = e) = 0. The configurational entropy is 
finite only for nontrivial solutions k > e. The tem- 
perature, Ta, where this happens for the first time is 
equal to the dynamic transition temperature of the sys- 
tem. The result for S c (T) also enables us to determine 
the Kauzmann temperature, Tk < 2a, where S c van- 
ishes. In Fig. 3, we show S C (T) for two different |e§| . The 
configurational entropy vanishes as Tk is approached: 
S c oc T — Tk- The regime (Ta — Tk)/Ta between the 
onset of slow dynamics and entropy crisis increases for 
decreasing strength of the phase separation. 



An important simplification of our approach resulted 
from the simple form, EqEI of the correlation function 
in replica space. All diagonal elements as well as all off 
diagonal elements are assumed to be identical. Whether 
this assumption is indeed stable can be addressed by eval- 
uating the eigenvalues of the stability matrix 



H, 



S 2 F 



a/3;7<5 



SG ab (q) SG^ 5 (q') ' 



(49) 



If there are negative eigenvalues of H our assumption 
for the replica structure is unstable. Following ReflHol 
we find that the lowest eigenvalue with respect to the 
replica indices is determined by the lowest eigenvalue of 
the matrijs^ 



K,* = <*(q - q') K ~ 2 (?) 



"0 



(50) 



5 $ 



In deriving this 



in momentum space, where vq 

with distinct a, /3, 7 and 5 

SG SG 

result we started from Eq[8] but used the fact that the 
functional <E> only depends on the momentum averaged 
correlation function, such that Vq becomes momentum 
independent. Ea l5UI is similar to the Schrodinger equa- 
tion in momentum space of a single particle with bare 
Hamiltonian K~ 2 (q) and local potential v$. The lowest 
eigenvalue, E, of this problem is given by 



d 3 
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(K- 2 (q)+E) 



(51) 



This equation can be analyzed if we find a way to calcu- 
late Vq which, due to EqO is determined by the first 
derivative of the self energy with respect to the cor- 



relation function 



5G-{5 



This derivative can be evalu- 



ated by following closely Baym and KadanoflS2. First 
we add to the local Hamiltonian an additional term 
— y] Q p Uupi'ai'P an( i analyze all correlation functions for 
finite U. At the end we will take the limit U — > and 
the correlation function has the simple structure Eg 1501 
For finite U , the correlation function is determined from 
the function Z [U] : 



G a (3 [U] = (4>a4>f3) 



6 log Z[U] 

SU a B 



(52) 



The self energy is a functional of only G [U] and not of 
U explicitly such that 



(53) 



If we furthermore introduce 



J a(3;j8 



SU S p 

{4>a<j>!34 > l4 l s) - (050/3) (54) 
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which can be evaluated explicitly once [7 = 0, one finds 
^2 { G al G ^ ~ ~$^T) L ^l3-n& = <W<W (55) 
which determines and thus vq. 

oG al 

Applying this approach to the Brazovskii model we 
find that the replica structure is marginally stable at 
the temperature Ta where the glassy state occurs for the 
first time. Below T\ the replica symmetric ansatz Eg 1501 
becomes unstable however it can be made stable if the 
replica index m does not approach 1 anymore but rather 
takes a value m = which defines the effective temper- 
ature T c g of the glass. This is a situation similar to one 
step replica symmetry breaking with break point given 
by mil A detailed discussion of this issue is presented in 
Ref|H 

Stability of the replica symmetric ansatz was only pos- 
sible because we consistently made the assumption of the 
dynamical mean field theory. Going beyond the local ap- 
proach, for example by using cluster DMFT techniques, 
enables one to study whether or not nonlocal phenomena 
change the replica structure of the theory. 



III. CONCLUSIONS 

Based on the recent development in the theory of 
strongly interacting electron systems we developed a dy- 
namic mean field theory for self generated glasses in a 
continuum field theory model. The key assumption of 
our approach, which applies to physical systems with 
short range higher order virial coefficients, is that glass 
formation is the consequence of predominantly local cor- 
relations. It is then possible to map the problem onto 
a purely local theory with same interaction and with 
a Gaussian part of the energy which is determined self 
consistently. The comparatively simple approach can 
easily be applied to multi component systems. Most 
importantly, recent developments in the cluster DMFT 
approach allow one to generalize this theory to include 
at least short range non-local effects. Owing to its lo- 
cal character a generalization of the present treatment 
to include spatially inhomogeneous states (correspond- 
ing with instantons) and thereby explicitly treat the mo- 
saic formation predicted in random first order transitions 
seems well within reach. 
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APPENDIX A: SELF-CONSISTENT 
SCREENING APPROXIMATION 

In this appendix we summarize the solution of the Bra- 
zovskii model within self-consistent screening approxima- 
tion of Refs l28ll3ll In Ref. a complete numerical so- 
lution of the self consistent screening approximation was 
developed which did not make any assumption with re- 
spect to the momentum dependence of the self energy. It 
was found that £ (q) does weakly depend on q for q < 2qo 
allowing to approximate S (q) ~ £(90) momentum in- 
dependent. This simplification was helpful to develop 
an analytical solution of the same problem, presented in 
Ref|3j Here we will show that the results of Refl3ll can 
easily be re-derived by performing a dynamic mean field 
theory where the local problem is solved within pertur- 
bation theory. 

It is useful to first summarize the problems which result 
in a perturbative treatment of the problem, something 
that can already be seen in the liquid state. Consider a 
given Feynman diagram for the self energy of a theory 
with interaction u<jf (r = 4 in our case). The number 
of internal lines is I n = i (nr — 2), whereas the number 
of internal momentum integrations of the original theory 
is L n = I n — (n — 1) = n (| — l). In a local theory we 
have for each running line a contribution q^G (i.e. we 
replace G (q) by its momentum averaged value). Thus, 
we obtain for a contribution to the self energy of order n 

~ „» (<$) Ln [qo 3 G\ In (Al) 

Using r = 4 and G of EqEOl it follows <5£(") ~ 

2 T L {'^) ■ Clearl y> in the limit of £ < x > where 

glassiness was found in Refs l28l3"ll the perturbation se- 
ries is a sum terms which are increasingly larger in magni- 
tude. e. strongly divergent. The self consistent screening 
approximation corresponds to the sum 

Sssc ~-X)(-l) B *sW oc -^y-V- (A2) 

n=l qae' 2 

Even though Esse = — <Zo e 1S small for small qae 2 , it is 
not clear at all whether this procedure leads to reliable 
results, making evident the need to use a method not 
based on the summation of diagrams. 

Using the replica approach for the local problem we 
obtain within the self consistent screening approximation 
for the local self energies: 

al^K = -2(D K K + D K F + D F K) 

ag£ F = -2D F F (A3) 
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with collective propagators 



propagator, Eq|23 we find for Hp = (k 2 — e 2 



D K 
D F 



2 



U F D Z K 



i + n F D K 



(A4) 



where Hk = K 2 + 2KF and II ^ = F 2 . This corresponds 
to the summation of bubble diagrams as summarized in 
detail in Ref|3l| and follows from the functional 



$[G1 



.-d 



TrlogZT 



(A5) 



with D a b = Dx^ab + Dp. Since we are solving a classi- 
cal (time independent) local problem we do not need to 
perform any momentum integrations in the evaluation of 
the diagrams. 

Using the expressions for the momentum averaged 



Hp = 



47T£ 



i np 



q ( l + e^- 1 ) 



(A6) 



where v — 
-1 



T?uq 2 al 



(2ir) 2 "^0"W Witl1 the choice fl = fl7T<?0 - 

l-7qQ L this agrees with the result of Ref l3ll analyzed in 
the strong coupling limit v^ 1 — + 0. In what follows we 
simply use ao — q^ 1 to be able to compare the results 
of the self consistent screening with the solution of the 
DMFT approach of this paper. In the limit v^ 1 — ► 
a solution k > e exists for the first time if e < £a — 
0.84. At the transition ka — 3.1ea- A detailed physical 
discussion of these results is given in Reff3jJ F° r weak 

/ x 2/5 

coupling it follows ea — 0.422 I — J with ka — 2.19£a- 
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